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Abstract. We investigate the anisotropic spin motive force in a-(BEDT-TTF) 2 l 3 , which 
is a multi-layered massless Dirac fermion system under pressure. Assuming the interlayer 
antiferromagnetic interaction and the interlayer anisotropic ferromagnetic interaction, we 
numerically examine the spin ordered state of the ground state using the steepest descent 
method. The anisotropic interaction leads to the anisotropic spin ordered state. We calculate 
the spin motive force produced by the anisotropic spin texture. The result quantitatively agrees 
with the experiment. 


1. Introduction 

An organic conductor (a-(BEDT-TTF)2l3 is a multi-layered massless Dirac fermion system, in 
which conduction layers of BEDT-TTF molecules and I3 anions stack alternatively. In each 
conduction layers, the massless Dirac fermion system is realized. The valence band and the 
conduction band contact at two inequivalent points in the Brillouin zone. The energy dispersion 
is linear in the vicinity of the contact points. The Fermi energy coincides with the energy of the 
contact point. Therefore, the system is called Dirac fermion system. 

Under the magnetic field S, the energy of the n-th Landau level of the Dirac fermion is 
described hy E — sgn{n)C^J\n\B with a constant C. The energy difference between the first 
and the zeroth Landau levels is ( 7 \/S. Therefore, in relatively high magnetic field, we may 
consider only the Landau level with n = 0 . 

(a-(BEDT-TTF)2l3 has unique spin ordered states [ 1 ]. In each layer, the quantum Hall 
ferromagnetic state is realized by the intralayer ferromagnetic interaction due to the exchange 
interaction. On the other hand, the interlayer ferrimagnetic state is possible by the interlayer 
antiferromagnetic interaction. 

(a-(BEDT-TTF)2l3 has anisotropic crystal structure and transport property. In figure [U the 
direction of the interlayer tunneling is shown. The tunneling is tilted in the direction of 6-axis. 
For the experimental data, the angle between the vertical direction and the interlayer tunneling 
direction is about 28 °. 

At low temperature and under magnetic field, anisotropic voltage is observed in (a-(BEDT- 
TTF)2l3 [2]. At high temperature the voltage disappears. The voltage increases with increasing 
the magnetic field. Therefore, one possible scenario is that the voltage is caused by spin motive 
force. The anisotropic crystal structure can lead to in the anisotropic voltage. The spin motive 
force produces the voltage. However, this is not a equilibrium situation. The spin motive force 
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Figure 1. The red arrow shows the interlayer tunneling direction. The BEDT-TTF molecules 
are described by the circles, a- and 6-axis is in the conducting plane, c-axis is perpendicular to 
it. (a) The interlayer tunneling direction in a-c plane. (6) The interlayer tunneling in b-c plane 
is tilted in the direction to the 6-axis. 


is transient due to the Gilbert damping. Charge accumulates at the edge of the system by the 
spin motive force and produces a electric field opposite to that created by spin motive force. 
The voltage will vanish due to the electric field created by the charge accumulation. Hereafter, 
for simplicity, we assume that the necessary time of the disappearance of the voltage is long and 
the voltage is constant in this time scale. 

In this paper, we show that anisotropic spin ordered states are realized in (a-(BEDT-TTF) 2 l 3 
because of the interlayer anisotropic interaction and anisotropic spin motive forces are produced 
under magnetic field. 
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Figure 2. The anisotropic voltage in (a-(BEDT-TTF) 2 l 3 . The voltage in the direction of a^b 
are defined by Va^Vi). 


2. Model 

We assume the square lattice for the conducting plane and the ferromagnetic interaction Jintra 
is assumed between nearest neighbor molecules. We define 6- and a-axis as x- and ^-axis, 
respectively. We take the direction of interlayer tunneling as z-axis. The interlayer interaction 
in the direction of tunneling is denoted by Jinter < 0. In this system, the lattice is distorted in 
x-z plane since the molecular arrangement is tilted in the direction of x-axis and there are two 
different next nearest neighbor sites in x-z plane. When the distance is different, the interaction 
should be different. For this reason, we assume a finite interaction J/^ter next 









nearest neighbor sites in x-z plane. On the other hand, for simplicity, we take a square lattice 
in x-z plane in this calculation and assume the magnetic field is parallel to the z-axis. 
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Figure 3. The model of oi-(BEDT- 
TTF) 2 l 3 . (a) Only Jintra exists in 
x-y plane. {h) Jinter is along 2 ;- 
axis and Jintra is along y-axis in y- 
z plane, (c) Jinter is along z-axis, 
Jintra is along x-axis and J/^ter i^ 
between one side of the next nearest 
neighbor sites in x-z plane. 


The Hamiltonian under the magnetic field is written as 


H = 


+ 


Jintra ^ ^ T ^i,jP-l,k T ^i—l,j,k T — 

Jinter ^ ^ {^iJ,k-\-l T • ^ij^k T *Anter ^ • ^ij^k 

ij^k hj^k 


(1) 


where ^ij^k is a spin at the lattice point (z, j, k). ^ is g factor and is the Bohr magneton. 

In this model, we investigate the spin ordered state of the ground state using the steepest 
descent method. At first, we put a spin with random direction on each lattice point. Next, we 
update spins with 
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Here, is the normalized localize moment. The resulting converged state is an 

approximate for the ground state. Jintra is proportional to the \fB [1], but, for simplicity, we 


























assume Jintra is constant. We numerically calculate with the number of sites, 12 x 12 x 12 and 
under the periodic boundary condition. 

The result with the parameters Jintra = — 5 K, Jinter = 4 K, J/^ter = — 2 K and B = b 
T is shown in figure 01 Here, we assume relatively large value for Jinter, which is the same 
order of magnitude as the interlayer hopping estimated in a related organic compound [3]. A 
ferromagnetic state is realized along y-axis. On the other hand, a periodic spin texture is created 
along x-axis due to the interlayer next nearest neighbor interaction. 



Figure 4. The resulting spin ordered state of (a-(BEDT-TTF) 2 l 3 with Jintra = ~5 K, Jinter = 4 
K, J/nter = “2 K and H = 5 T. (a) The projection of spins in a layer to x-y plane. A periodic 
spin ordered state is realized along x-axis, although a ferromagnetic structure is along y-axis. 
{h) The projection of spins in a layer to x-z plane. 


Assuming this spin ordered state, we calculate the spin motive force [4] in the strong coupling 
limit. In this limit, the spins are polarized and there are no conduction electron with minority 
spin. Therefore, the spin motive force is identical with the voltage. We make indices j, k 
implicit below since we consider spin motive force in uniaxial direction. The localized moment 
is represented as 


= (cos cos 0^, cos sin 0^, sin 0^), (3) 

where 6i is the angle between the direction of the spin and the z-axis. (pi is the angle between 
the direction of the spin projected in the x-y plane and the positive direction of the x-axis. The 
electric field due to the localized moment is written as 

Ei = A sin 0. [{dtOi) (V0,) - {dtct>i) (V0i)], (4) 

where e is the charge of the electron [5]. Here, the condition dtOi = 0 and dffpi = g/j^BBlh are 
always satisfied since the magnetic field is along z-axis. The electric field is rewritten as 

= (5) 

In the numerical calculation with the periodic boundary condition, the electric field is canceled 
if we take the average over the spatial period and so the spin motive force does not appear. 
However, in the real system, there is the edge of the sample and a finite spin motive force is 
possible. On the other hand, spin motive force along ?/-axis does not appear at all. Therefore, 




the spin motive force in this model has strong anisotropy. In order to estimate the spin motive 
force in the real system, we calculate the spin motive force with taking a part of the periodic 
spin structure. 

The numerical result is shown in figure [5l The spin motive force along x-axis is on the order of 
0.01 mV to 0.1 mV which is about the same order compared with the experimental result. The 
most of electric field is canceled and only the spin motive force in the vicinity of the edge remains. 
The experimental result shows that increases with increasing the magnetic field below 7 T. 
In order to investigate the magnetic field dependence of we define the maximum value of the 
Vb as Vnax- Although Vnax docs not always correspond to the experimental observed voltage, it 
is able to qualitatively evaluate the magnetic field dependence of Vb. The nurerical result of the 
magnetic field dependence of Vnax is shown in figure El Vnax increases with increasing magnetic 
field below 4 T. However, above 5 T, Vnax decreases. Under high magnetic field, the Zeeman 
energy is larger than the intralayer and the interlayer interactions, and so {V0^} are small 
since we assume the interactions are constant. Therefore, the electric fields due to the localized 
moments vanish under high magnetic field. If we consider the magnetic field dependence of the 
interactions, we would evaluate the anisotropic spin motive force more precisely. 
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Figure 5. Spin motive forces along x-axis in (a-(BEDT-TTF) 2 l 3 . rix is the number of spins 
along x-axis from the edge of the system. 

In a real system, the edge shape is different for a different layer. So, the naive expectation 
is that if we take the average over those edges, the spin motive force is canceled out. However, 
since the anisotropic spin motive force is generated by the three-dimensional structure, spins in 
regions without interlayer interaction should not contribute to the spin motive force. Therefore, 
it is justified to consider spin motive force for layers with the same edge shape. 

3. conclusion 

We have calculated the spin motive force in (a-(BEDT-TTF) 2 l 3 . The anisotropic interlayer 
interaction exists reflecting the anisotropic crystal structure. The anisotropic interaction leads 
to the unidirectional periodic spin structure. The anisotropic spin motive force is created from 
the anisotropic spin structure under magnetic field. 
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Figure 6. The magnetic field dependence of the maximum values of V 5 . Idnax increases with 
increasing magnetic field below 4 T and decreases above 5 T. 


Acknowledgments 

We would like to thank N. Tajima for discussions and sending us the experimental result. 
This work was financially supported in part by a Grant-in-Aid for Scientific Research (A) on 
Dirac Electrons in Solids (No. 24244053) and a Grant-in-Aid for Scientific Research (B) (No. 
25287089) and (C) (No. 24540370) from the Ministry of Education, Culture, Sports, Science 
and Technology, Japan. 

References 

[1] Kubo K and Morinari T 2014 J. Phys. Soc. Jpn. 83 033702 

[2] Tajima N Private communication 

[3] Jindo R, Sugawara S, Tajima N, Yamamoto M, Kato R, Nisino Y and Kajita K 2013 Phys. Rev. B 88 075315 

[4] Tatara G, Kohno H and Shibata J 2008 Physics Reports 468 213 

[5] Tserkovnyak Y and Mecklenburg M 2008 Phys. Rev. B 77 134407 




